arXiv: 1506.07045v5 [hep-ph] 7 Apr 2017 


Production of or and its excited states via t-qnark or t-quark decays 

Qi-Li Liao ^0 Jia-Wei Zhang ^0 Guo-Ya Xie Yan Yu Ya Deng and Guang-Chuan Wang ^ 

^ College Mobile Telecommunications Chongqing University of 
Posts and Telecom, Chongqing 401520, People’s Republic of China 
^ Department of Physics, Chongqing University of Science and Technology, Chongqing 401331, People’s Republic of China 
® Chongqing University of Posts and Telecom, Chongqing 400065, People’s Republic of China 

(Dated: April 10, 2017) 

In this work we evaluate the masses of the |(6s)[n-]) or |(fes)[n]) quarkonium (73° or meson) 
under the B.T. potential, and the values of the Schrodinger radial wave function at the origin of the 
|(?)s)[n]) or |(&s)[n]) quarkonium within the five potential models. Then we investigate a systematic 
study on the production of the |(fes)[n]) or |(6s)[n]) quarkonium via top quark or antitop quark 
decays in the color-singlet QCD factorization formula (CSQCDFF), i.e., the two 5'-wave states, 

|(6s)[1^5'o]) (or |(6s)[l^S'o])) and |(fos)[l^S'i]) (or |(6s)[l®S'i])), and its four P-wave excited states, 
i(64[l^Pi]) (or |(6s)[DPi])) and |(te)[PPj]} (or |(Bs)[PPjj}) (with J = [0,1,2]). For deriving 
compact analytical results for complex processes, the “improved trace technology” is adopted to 
deal with the decay channels at the amplitudes. Moreover, various differential distributions and 
uncertainties of the concerned processes are analyzed carefully. By adding the uncertainties caused 
by the b and s-quark masses in quadrature, we obtain r(t —>■ |(6s)[n]} + W^s) = Id.lOi'ig'jo MeV. 

At the LHC with the luminosity C oc 10®^and the center-of-mass energy '/S = 14 TeV, 
sizable |(&s)[n-]) or |(bs)[n]) meson events can be produced through t-quark or t-quark decays; i.e., 
about 1.3 X 10® P® or P® events per year can be obtained. 

PACS numbers: 12.38.Bx, 14. 65.Ha, 14.80.Bn, 14.40.Pq 


I. INTRODUCTION 

The GDF and DO Gollaborations at the Tevatron have successfully collected numerous data for the P® and P® 
meson The LHG shall also provide a good platform to study the properties of the P® and P® meson, e.g., recent 
results and review papers on the P® and P® meson measurements and searches by the LHCb, GMS, and ATLAS 
Gollaborations at the LHG can be found in Refs. 1M3. 

The P® and P® meson form an interesting experiment and theory for the study of the quantum chromodynamics 
(QGD) and new physics (NP), which can be used to study those interesting topics as the QCD model building, the 
electroweak symmetry breaking mechanism, charge-parity (CP) violation, new physics beyond the Standard Model 
(SM) and etc. |lll - [ll| . Taking into account the prospects of the P® or P® meson production and decay at the 
TEVATRON and at the running LHC, the future numerous data require more accurate theoretical predictions on its 
hadronic production and indirect production through weak-decay processes. 

For the P® or P® meson direct production, various investigations has been carried out studied in detail in Refs. di¬ 
ll^. The P® or P® meson indirect production through the Z® boson and high excited T(nS) meson decays has 
been measured and discussed in Refs. (lOl - l^ . Alth oug h in literature the direct hadronic production of the P® or 
P® meson has been thoroughly studied in Refs. [l3, Il8l | and CDF discovered the meson which just comes from the 
hadroproduction. As a compensation to understanding the production mechanisms, and also testing the perturbative 
quantum chromodynamics (pQCD) and the QCD factorization formula. It is quite interesting that to study the 
indirectly production of the P® or P® through top quark or antitop quark decays, especially on considering the fact 
that numerous t-quark or t-quark may be produced at the LHC. At the LHC running at the center-of-mass energy 
'/S = 14 TeV and luminosity raising up to £ oc [H, [2^, about 10® ti quark per year will be produced. 

More t-quark or t-quark rare decays can be adopted for precise studies at the LHC. This paper is devoted to study the 
indirect production of the P® or P® meson via t-quark or t-quark decays. The theoretical studies of the bound states is 
based on the factorization formula, which includes the color-singlet QCD factorization formula (CSQCDFF) (2^, the 
nonrelativistic quantum chromodynamics (NRQCD) [2^, the fcj^-factorization [2^ and etc. The CSQCDFF is a useful 
method for studying on the production and decays of the quarkinum. According to the CSQCDFF, the processes of 
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the production and decays of the quarkinum can be factorized as two parts of the perturbatively calculable short- 
distance coefficients and the nonperturbative long-distance factors. In present paper it is investigated a systematic 
study of the indirect production of the |(6s)[n]) or |(6s)[n]) quarkonium via t-quark or t-quark decays within the 
CSQCDFF. To derive the analytical expression of the pQCD calculable conventional squared amplitudes becomes too 
complex and lengthy for more (massive) particles in the final states and for the excited states to be generated for 
the emergence of massive-fermion lines in the Feynman diagrams, especially to derive the squared amplitudes of the 
P-wave states. In Refs, [sol - ls^ . the “improved trace technology” is suggested and developed to solve the problem. 
Here we adopt the “improved trace technology” to derive the analytical expression for the above-mentioned processes. 
Without confusing and for simplifying the statements, later on we will not distinguish the and meson unless 
necessary. Because of the production of the |(^^s)[n]) or |(&s)[n]) meson through t-quark or t-quark decays channels, 
i.e., t |(^^s)[n]) -I- W^s and i |(^s)[n]) -I- W~s, and t —>■ |(&s)[n]) -I- VF+6 and i —|(6s)[n]) -I- W~b, are symmetric 
in the interchange from particle to anti-particle. 

The rest is organized as follows. In Sec. II, we introduce the calculation techniques for the above-mentioned t-quark 
and t-quark decays to the |(&s)[n]) quarkonium. In order to calculate the production of the |(6s)[n]) quarkonium and 
its excited states, we evaluate the masses of the |(6s)[n]) quarkonium with the energy eigenvalue of the nonrelativistic 
Schrbdinger equation. In Sec.Ill, we calculate and tabulate the masses of the |(l)s)[n.]) quarkonium and its values of 
the Schrbdinger radial wave functions at the origin. Then, we investigate a systematic study on the production of 
|(6s)[n]) through t-quark and t-quark decays channels, i.e., t —>■ |(6s)[n]) -I- kF+s and t — >• \{bs)[n])+W~b, where |(6s)[n]) 
stands for |(6s)[l^S'i]), |(6s)[l^S'i]), |(6s)[l^Pi]), and |(6s)[l^Pj]) (with J = [0,1, 2]). We make some discussions on the 
various differential distributions and the uncertainties of the decays widths by the masses of the |(6s)[n]) quarkonium 
and the correspondingly values of the Schrbdinger radial wave functions at the origin. The fourth section is reserved 
for a summary. The idea of “improved trace technology” is put in the appendix finally. 


II. CALCULATION TECHNIQUES AND MATRIX ELEMENT {O^(n)) 

A. Calculation techniques 

For the production of the |(6s)[n]) meson through t-quark and t-quark decays channels, t{k) —>■ |(l's)[n])(g 3 ) -I- 
W^{q 2 ) + s{qi) and i{k) —>• |(l's)[n.])(g 3 ) -I- W~{q 2 ) -f &(gi), where k and qi {i = 1,2,3) are the momenta of the 
corresponding particles, according to the CSQCDFF, its total decay widths dF can be factorized as 

dr = ^df(t^ |(6s)[n])-f VF+s)(0^(n)), 

n 

dr = ^df(t ^ |(6s)[n])-f VF-6)(Ci^(n)). (1) 

n 

where the nonperturbative matrix element {O^ in)) describes the hadronization of a bs pair into the observable quark 
state H and is proportional to the transition probability of the perturbative state bs into the bound state |(ds)[n]) 
quarkonium. The matrix elements {O^(n)) for the color-singlet components can be directly related to the Schrbdinger 
wave functions at the origin for the S'-wave states and the first derivative of the wave functions at the origin for the 
P-wave states, which can be computed via the potential mod els J37H42| . The nonperturbative matrix element {O^(n)) 
can also be obtained via potential NRQCD (pNRQCD) l44lj aird/or lattice QCD [d^. 

The hard-scattering amplitudes for specified processes can be dealt with 

t —|(&s)[n]) -I- 

i ^ \{bs)[n\)+W~b. (2) 

The Feynman diagrams of the two processes, t{k) —>• |(^^s)[n])(g 3 ) -I- W^{q 2 ) -f s{qi) and t{k) —|(&s)W)((Z 3 ) + 
W~{q 2 ) + b{qi), are presented in Figs. [1] and [2J where the intermediate gluon should be hard enough to produce a ss 
pair or bb pair, so the amplitudes are pQCD calculable. 

These amplitudes of the above-mentioned channels can be generally expressed as 

m 

iM = CsUsiiqi) ^ AnUs'jik), (3) 

n—1 

where m stands for the number of Feynman diagrams, s and s' are spin states, and i and j are color indices for the 
outing quark and the initial top quark or antitop quark, respectively. The overall factor Cs = {2gg‘fVtq/3'/2)5ij stands 
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FIG. 1 : (color online). Feynman diagrams for the process t{k) |(6s)M)(93) + W'^{q2) + s{q\), where |(fos)[n]) stands for 
|(fos)[^ 5 o]), |( 6 s)p 5 'i]), |(6s)[^Pi]), and |(6s)pPj]} (with J = [ 0 , 1 , 2 ]), respectively. 



FIG. 2 : (color online). Feynman diagrams for the process t{k) —>■ |(fos)[n])(g3) + W (52) + b{qi), where |(6s)[n]) stands for 
|( 6 s)[^ 5 'o]), |(&s)[®Si]), |(6s)[^Pi]), and |(&s)pPj]) (with J = [ 0 , 1 , 2 ]), respectively. 


for the specified quarkonium in the color-singlet, here Vtq is the Cabibbo-Kobayashi-Maskawa (CKM) matrix element, 
where Vtq = Vtb stands for the processes t —>• |(6s)[n]) -|- FF+s, and Vtq = Vts for i —>• |(6s)[n]) -I- W~h. While for 
production the (^-Si) states of the (6s)[n]) quarkonium via the decays channel t{k) —>• |(6s)[n])(q3)-|-FF+(g2) + s(9i), 
the amplitudes An can be written as 


Ai — 


A2 — 


(91+ 932)2 ^“(<11+< 73)27 ) 


7 n 7 -^^-^/( 92 )(l- 7 ^) 


53 v/<2 + ^ 3 l)+""‘ 


' (< 7 i + 932)2 

For the ^Pi-wave states. An can be written as 


(92 + 931)2 - mt 


q =0 


q =0 


^f=o+=i = efiqs) 


dQfi 


^S'=0,Z/=1 (1 / \ 

A = q (93)1— 


A 

' dq^ 


(91 + 932 ) 2 ^“ ( 91 + 93 ) 27 ) 
^ ng3(g) _u, VI „5 /^2 + ^3l) + ""* 


and the ^Pj-wave states (J = 0 , 1 , 2 ) 


. 5 = 1 ,L =1 J f \ ^ 

dq, 


'“(gi -h 932)2'"^'"'^^" ' ^(g2-h 931)2 -TO' 

As A A + 


q =0 


q =0 


7“?—r^^7c<-pV—ry-- tA92)(1 - 7*^) 

( 91 + 932)2 (91 + 93 ) 2--"2 


. 5 = 1 ,L =1 

^2 




= A (qs) 


la 




(91 + 932)2' 


(92 + 93i)2 _ ^2 


J 9=0 


g=0 


( 4 ) 

( 5 ) 


( 6 ) 

( 7 ) 


( 8 ) 

( 9 ) 


Here q stands for the relative momentum between the two constituent quarks in the |(&s)[n]) quarkonium. 931 and 
932 are the momenta of the two constituent quarks, i.e., 


mb ms 

931 = - 93 + 9 , 932 = - 93 - 9 - 

Tfibs T^bs 


( 10 ) 


Where mbs = mb + ms is implicitly adopted to ensure the gauge invariance of the hard scattering amplitude. £(92) is 
the polarization vector of W'^ boson. £5(93) and £((93) are the polarization vectors relating to the spin and the orbit 
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angular momentum of the |( 6 s)[n]) quarkonium, S'^^{q 3 ) is the polarization tensor for the spin triplet P-wave states 
(with J = [0,1, 2]). The covariant form of the projectors can be conveniently written as 

+ “‘I <8) (“) 

and 

+ ( 12 ) 

Here Ic stands for the unit color matrix with Nc = 3 for the QCD. 

After substituting these projectors into the above amplitudes and doing the possible simplifications, the amplitudes 
then can be squared, summed over the freedoms in the final state and averaged over the ones in the initial state. The 
selection of the appropriate total angular momentum quantum number is done by performing the proper polarization 
sum. If defining 




<l3fj.Q3v 

M 2 


(13) 


the sum over polarization for a spin triplet 5'-state or a spin singlet P-state is given by 


(14) 


where Jz = Sz or Iz respectively. In the case of ^Pj states, the sum over the polarization is given by [46 


p.( 0 )p( 0 )* 

Ai)p(i)* _ 






— 2 H^i/zH^/^y) 


Vf(2)p(2)* ^ 




(15) 

(16) 

(17) 


The amplitudes An of the t —>■ |(6s)[n]) + W~b in the formulas are similar in the above-mentioned decays channels, 
so do not listed here to shorten the paper. 

Finally, the short-distance decay widths are given by 


df(t^|( 6 s-)[n]) + lF+s) = 


2 fc 0 ' 


dr{t^\{bs)[n]) + W-b) = _^|Mpd$ 3 . 


(18) 


where ^ means that one needs to average over the spin and color states of the initial particle and to sum over the 
color and spin of all the final particles. 

In the t-quark or t-quark rest frame, the three-particle phase space can be expressed as 


1^'= - E ?/j n (4%- 

The 1 —>■ 3 phase space with massive quark/antiqark in the final state have been calculated to simplify procedures 
in Refs. [^ . 1^ . To shorten the paper, we shall not present it here and the interested reader may turn to these 
two references for the detailed technology. We can not only derive the whole decay widths but also obtain the 
corresponding differential decay widths that are helpful for experimental studies with the help of the formulas listed 
in Refs. [13, such as dT/dsi, c?r/ds 2 , dT/dcos9i2, and c?r/dcos 0 i 3 , where si = {qi +q 2 )^, S 2 = (91 + 93 )^, ^12 
stands for the angle between 91 and 92 , and 613 for 91 and 93 . 

For the present considered the squared amplitude of top quark semi-exclusive decay process, t{k) —>■ ( 6 s)[n]( 93 ) -I- 
hF~'"( 92 ) + s( 9 i), there are two massive spinors of top quark and s-quark. So to achieve the analytical expression for the 
squared amplitude becomes too complex and lengthy in the conventional way, especially for P-states of the |(5s)[n]) 








5 


meson. For solved the problem, we adopt ‘new trace amplitude approach’ to deal with the hard scattering amplitude 
( 151 ) for such complicated processes at the amplitude level. In a difference from the helicity amplitude approach [13- 
l50j, only the coefficients of the basic Lorentz structures are numerical at the amplitude level. However, by using the 
“improved trace technology” in Refs, one can sequentially obtain the squared amplitudes, and the numerical 

efhciency can also be greatly improved. Detailed processes of the approach can be found in Refs. and here 

for self-consistency, we shall present its main idea in appendix. 

Finally, the partial decay widths over si and S 2 can be expressed as 


dr 

dsids2 


2567r3m? 




( 20 ) 


where rrit is the mass of the top quark. 

According to the CSQCDFF [^, the color-singlet nonperturbative matrix element {O^ (n)) can be related to the 
Schrodinger wave function (0) at the origin and the first derivative (®) ^ P-wave states of the | {bs) [n]) 

quarkonium. 

(O^inS)) |V^|(fc,)[„5]>(0)P, 

{O^inP)) (21) 

As the spin-splitting effects are small, we will not distinguish the difference between the wave function parameters 
for the spin-singlet |(6s)[l^S'i]) and spin-triplet |(6s)[l^S'i]) states, also not discriminate the difference among the wave 
function parameters for |(6s)[l^Pi]) and |(6s)[l^Pj]) (with J = [0,1,2]) in following calculation. 


B. Mass of the |(&s)[n]) quarkonium and matrix element {O^(n)) 

The mass spectrum of the |(6s)[n]) bound states has the form 

M„(|(6s)[n])) = mb + rris + £’„(mt,, ms,V{fj). (22) 

where mt, and are the current quark masses in the MS scheme, and En{mb, mg, H(r )) is the energy eigenvalue of 
the nonrelativistic Schrodinger equation with a flavor-dependent potential models V (r) [371 ] . 

According to the CSQCDFF [^, the color-singlet nonperturbative matrix element {0^{n)) can be related to the 
wave function at the origin. In the rest frame of the |(&s)[n]) quarkonium, it is can be separated the Schrodinger wave 
function into radial and angular pieces 'I'nZm(fO = Rni{r)Yim{d,‘p), where n is the principal quantum number, and I 
and m are the orbital angular momentum quantum number and its projection. Rni{r) and Yim{0,^) are the radial 
wave function and the spherical harmonic function accordingly. The wave function at the origin 4'|f,s)[raS])(0) the 
first derivative (0) are related to the radial wave function at the origin P|(bs)[riS])(0) and the first derivative 

^I(&i)[nP]>(0)^ accordingly. 


'I'|(&s)[nS]>(0) = Vl/47ri?|(;,s)[„5])(0), 

'^'Kbs)[nP]>(0) = VV^^\(bs)[nP])i^)- (23) 

III. NUMERICAL RESULTS AND DISCUSSIONS 
A. Input parameters 

The values of the wave function at the origin for the |(6s)[n]) quarkonium are related to the number of flavor 
quark n/, the Regge slope a', the 1-loop beta-function coefficient /3o and the 2-loop beta-function coefficient /3i, 
the renormalization scale parameters and so forth, where MS is the modihed minimal subtraction scheme. In 

calculating the wave function at the origin of the |(6sQ[n]) quarkonium within the five potentials Him, we adopt the 
Regge slope a' = 1.067GeU“^, and the scale parameters A-^ as A!!lr^=0.386 GeV, Aff£r^=0.332 GeV, Aff£r^=0.231 
GeV [^. The current quark masses, i.e., mb = 4.19 GeV and mg = 0.101 GeV, are chosen as derived in Ref. 
for calculating masses of the |(6s)[n]) quarkonium. The masses of the |(6s)[n]) quarkonium {B^ meson) under the 
Buchmialler and Tye potential (B.T. potential) [13,113, [H, 13 are presented in Table HI that can be derived through 
solving the energy eigenvalue of the Schrodinger equation [5l| . The values of the constituent quark masses mb and 
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TABLE I: Mass of the |(6s)[n.]) quarkonium under the B.T. potential (^. 1^. 



n = 1 

n = 2 

n = 3 

n = 4 

n = 5 

S states (GeV) 

5.37 

6.08 

6.70 

7.18 

7.60 

P states (GeV) 

5.83 

6.47 

7.02 

7.53 


D states (GeV) 

6.15 

6.75 

7.29 



F states (GeV) 

6.50 

7.06 





rus for the |(fe)[n]) quarkonium and the accordingly quantities |i?|(f,s)[„s]>( 0 )P> l^|(f,s)[nP])( 0 )p! and |[„£,])( 
are presented in Table mi for the five potential models [37|. The five potentials are the B.T. potential, the John L. 


Richardson potential (J. potential) 


the K. Igi and S. Ono potential(I.O. potential) [s^, the Yu-Qi 

Chen and Yu-Ping Kuang potential(C.K. potential), and Coulomb-plus-linear potential (the so-called Cor. potential) 
[13) HI) HI) 1131 , respectively. During the following calculation, we adopt the values of wave functions at the origin 
and the constituent quark masses nib and nig for the |(&s)[n.]) quarkonium under the B.T. potential in Table |TT] (with 
n/ = 3) as the central values. 

The input other parameters are adopted as the following values [12, [11|: tou/= 80.399 GeV, mt = 172.0 GeV, 
|Vtt,|=0.88, |V4 s|=0.0406, Gf = 1.16639 x 10“^ GeV~^. We set the renormalization scale to be niBs of the \{bs)) 
quarkonium for leading-order as running , which leads to Q;s(niSs)=0.26. Furthermore, we adopt the constituent 
quark mass nig = 0.50 GeV and nib = 4.87 GeV for the |( 6 s)[lS']) states, and mg = 0.69 GeV and mb = 5.14 GeV 
for the I( 6 s) [IP]) states. As the masses of s-quark is greater than Ajjg, so the above-mentioned processes is pQGD 
calculable. To ensure the gauge invariance of the hard amplitude, the |( 6 s)[n]) quarkonium mass M is set to mb + mg, 
where mb and mg stand for constituent quark mass. 


B. production of the |(6s)[n]) quarkonium via t-quark and t-quark decays 


As a reference for branching fractions, we calculate the decay widths for the basic processes t ^ b + W~^ and 
i ^ s + W~. Their decay widths can be written as 


r 


\Vt,\^GF\p\ 

2 y/ 2 Trmt 


+ |p|2 -k2|p|2( ^Jml + jpP + \j'm\ -k jpP ) 


(24) 


where Vtq = Vtb and mq = mb are for processes t —>■ 6 J- W'^, and Vtq = Vts and mq = mg ioi t ^ s + W . p stands 
for the relative momentum between the final two particles in the rest frame of the top quark or antitop quark. 


1^1 ^ V{m'i - {mq - mwY){ml - {mq + mwY) ^ 25 ) 

2 mt 

Then, we obtain T{t b + W^) = 1.131 GeV and r{t s + W~) = 2.416 MeV. 

The decay widths and branching fractions for the |( 6 s)[n]) quarkonium states through the two decay channels, 
t —>• |( 6 s)[n]) -I- IF+s and t —>• |( 6 s)[n]) J- W~b, are listed in Tables IIIII and HVl within the B.T. potential. 

To show the relative importance among different states more clearly, we present the differential distributions dT/dsi, 
dT/ds 2 , dT/dcos9i2, and dT/dcosOis for the above-mentioned channels in Figs. [3] and S] Moreover, we define a ratio 


_ (ir/dsi(|(6s)[n])) 

" dr/ds,(|(6s)[D5o]))’ 


(26) 


where i = 1, 2, and n stands for l^Pi, and l^Pj (with J = [0,1, 2]), respectively. The curves are presented in Fig. [5] 
These figures show explicitly that the excited states, i.e., |( 6 s)[l^Pi]) and |( 6 s)[l^Pj]) (with J = [0, 1, 2]), can provide 
sizable contributions in comparison to the lower state |( 6 s)[l^S'o]) or |( 6 s)[l^S'i]) in almost the entire kinematical 
region. 

For the production of |( 6 s)[n]) through channel t —>• |( 6 s)[n]) + sW~^, its total decay width for all the P-wave states 
is 793.6 KeV, which is about 18% (9.2%) of the decay width of the |(6s)[1^S'q]) (|(6s)[1^S'q])) meson production. For 
the |( 6 s)[n]) production via i —>■ |( 6 s)[n]) J- bW~, its total decay width for all the P-wave states is 2.092 eV, which is 
about 29% (23%) of that of |(6s)[l^S'o]) (|(6s)[l^S'o])). But though t —^ |( 6 s)[n]) -I- bW~ is the GKM matrix element 
suppressed to t —>■ |( 6 s)[n]) + sVF+ by jits^ 0.2%, it is enhanced by the phase space, and it is easier to a 
gluon hadronic generated a ss pair than a 66 pair. So as a combined result, the decay width of t —>■ |( 6 s)[n]) + bW~ is 
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FIG. 3: (color online). Differential decay widths dT/ds\ and dV/ds 2 for t —>■ |(bs)[n]} + sW^ {nf = 3), where the dashed line, 
the dash-dotted line, the solid line, the crosses line, the diamond line, and the dotted line are for |(tis)[1^5o]), |(tis)[l®5'i]), 
|(fes)[l^Pi]), |( 6 s)[l^Pol), |(fos)[l®Pi]), and |( 6 s)[l®P 2 ]), respectively. 




FIG. 4: (color online). Differential decay widths dV/dcos 6\2 and dV/dcos 6 \i for t —>■ |( 6 s)[n]) -|- sW^ {uf = 3), where the 
dashed line, the dash-dotted line, the solid line, the crosses line, the diamond line, and the dotted line are for |(hs)[1^5'o]), 
|(bs)[l®5'i]), |( 6 s)[l^Pi]), |( 6 s)[l®Po]}, |( 6 s)[l®Pi]}, and |( 6 s)[l®P 2 ]), respectively. 




FIG. 5: (color online). The ratios P[l] and P[2] versus si and S 2 for the channel t — >■ |(fos)[n]) -I- sW^ {uf = 3). Here the dash- 
dotted line, the dotted line, the solid line, and the dashed line are for |(bs)[l^Pi]), |( 6 s)[l®Po]), |(&s)[l®Pi]), and |( 6 s)[l^P 2 ]), 
respectively. 
















TABLE II: Radial wave functions at the origin and related quantities for the |(fos)[n]) quarkonium. 


|(fos)[n]) 

Mass and potential 

n = 1 

n = 2 

n = 3 

n = 4 

n = 5 


mslrrih (GeV) 

0.50/4.87 

0.73/5.35 

0.98/5.72 

1.22/5.96 

1.46/6.14 


B.T.(n/=3) [3^ 

0.6516 

0.6737 

0.6889 

0.7025 

0.7808 


B.T.(n/=4) [3^ 

1.098 

0.8043 

0.7569 

0.8035 

0.9311 


B.T.(n/=5) [3^ 

0.7721 

0.7884 

0.7484 

0.7058 

0.8373 

S states 

J. (n/=3) 

0.5505 

0.6639 

0.8248 

0.9730 

1.116 


J. (n/=4) 

0.4814 

0.5667 

0.6978 

0.8185 

0.9346 


J. {nf=5) 

0.3169 

0.3559 

0.4329 

0.5053 

0.5760 

|-R|[nS]>(0)|^ 

1.0. (n/=3) Mi 

0.2801 

0.3497 

0.4407 

0.5241 

0.6041 

(GeV^) 

1.0. (n/=4) M 

0.2938 

0.3616 

0.4536 

0.5382 

0.6193 


1.0. (n/=5) 

0.2996 

0.3656 

0.4572 

0.5417 

0.6230 


C.K.(n/=3) M 

0.3444 

0.3911 

0.4745 

0.5521 

0.6276 


C.K.(n/=4) [42] 

0.3703 

0.4120 

0.4958 

0.5734 

0.6486 


C.K.(n/=5) M 

0.4042 

0.4383 

0.5221 

0.5997 

0.6741 


Cor. [38] 

0.3986 

0.5109 

0.6608 

0.8016 

0.9396 


rUs/mj, (GeV) 

0.69/5.14 

0.99/5.48 

1.21/5.81 

1.40/6.13 



B.T.(n/=3) [3^ 

3.029 

17.10 

26.21 

35.29 



B.T.(n/=4) [3^ 

13.32 

14.46 

37.13 

54.23 



B.T.(n/=5) [3^ 

8.568 

13.04 

25.24 

39.59 


P states 

J. (n/=3) M 

7.889 

22.23 

38.81 

57.60 



J. (n/=4) M 

6.036 

16.92 

29.39 

43.46 



J. (n/=5) [^ 

2.747 

7.324 

12.58 

20.40 


|-^|[raP]) (0)P 

1.0. (n/=3) M 

2.546 

7.377 

13.11 

19.69 


(10-^ GeC®) 

1.0. (n/=4) [^ 

2.692 

7.764 

13.75 

20.59 



1.0. (n/=5) [^ 

2.785 

7.974 

14.06 

20.99 



C.K.(n/=3) M 

3.362 

9.396 

16.22 

23.83 



C.K.(n/=4) M 

3.632 

10.13 

17.44 

25.56 



C.K.(n/=5) [42] 

3.993 

11.09 

19.03 

27.82 



Cor. M 

3.887 

11.59 

20.97 

31.93 



mafmb (GeV) 

0.75/5.40 

0.98/5.77 

1.18/6.11 




B.T.(n/=3) [3^ 

0.5089 

4.272 

12.77 




B.T.(n/=4) [3^ 

1.332 

12.93 

29.50 




B.T.(n/=5) [^ 

1.672 

7.200 

18.93 



D states 

J. (n/=3) M 

3.628 

15.17 

37.38 




J. (n/=4) [^ 

2.384 

9.940 

24.41 




J. (n/=5) [^ 

0.5959 

2.344 

5.443 



l%„n]>(0)P 

1.0. (n/=3) M 

0.7639 

3.257 

8.145 



(10"= GeV’^) 

1.0. (n/=4) M 

0.8089 

3.443 

8.591 




1.0. (n/=5) M 

0.8430 

3.572 

8.883 




C.K.(n/=3) [42] 

1.023 

4.307 

10.63 




C.K.(n/=4) M 

1.101 

4.637 

11.44 




C.K.(n/=5) M 

1.206 

5.081 

12.53 




Cor. M 

1.295 

5.602 

14.16 




very smaller than that of t —>■ |(6s)[n]) + sW~^ by only about 1 x 10 Thus, we will do not discuss the uncertainties 
of the |(6s)[n]) meson production through i —>■ |(6s)[n]) + bW~ in the next subsection. 

Running at the center-of-mass energy y/S = 14 TeV and with luminosity at the LHC, one may 

expect to produce about 10® tt-quark per year [2^ [2^. We can estimate the events number of the \{bs)) quarkonium 
production through t-quark and t-quark decays, i.e., 4.1 x 10^ |(6s)[1^5'o]) or |(6s)[1^5'o]), 7.8 x 10® |(6s)[l®5'i]) or 
|(6s)[l®S'i]), 7.0 X 10^ |(6s)[lP]) or |(&s)[lP]) events can be obtained per year. It might be possible to hnd P° or 
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TABLE III: Decay widths and branching fractions for the production of the |(6s)[n]) quarkonium through t-quark semiexclusive 
decays within the B.T. potential (n/ =3). 


t — >• |(6s)[n]) -1- W^s 

r(KeV) 

r(t-»|(6s)[nl)+W+s) 

r(t-»w++6) 

|(6s)[iSo])-fIF+s 

4528 

4.0 X 10“® 

|(6s)pSi])-kIF+s 

8614 

7.6 X 10“® 

|(6s)[iPi]) +IF+S 

293.3 

2.6 X lO”"* 

|(&s)pPo]) +IF+S 

119.7 

1.1 X lO”"* 

|(&s)pPi]) +IF+S 

142.5 

1.3 X lO”"' 

t^\ibs)[^P2]) +W+S 

238.1 

2.1 X lO""* 


TABLE IV: Decay widths and branching fractions for the production of the |(&s)[n]) quarkonium through t-quark semiexclusive 
decays within the B.T. potential (n/ = 3). 


t —>• |(fos)[n]} + W b 

r(eV) 

r(t— >-|(6s)[n]) + VK b) 
r(t^w-+s) 

|(te)[iSol}-tlF-b 

9.056 

3.7 X 10“® 

|(6s)p5il)-tlF-b 

7.095 

2.9 X 10“® 

t^\{bs)[^Pi]) + W-b 

0.560 

2.3 X 10“^ 

t^\{bs)[^Po]) + W-b 

0.915 

2.4 X 10“'^ 

i^\{bs)fPi]) + W-b 

0.531 

2.2 X 10“^ 

t^\{bs)fP2]) + W-b 

0.086 

3.6 X 10“® 


through top quark or antitop quark decays, since one may identify these particles through their cascade decay 
channels, e.g., —>■ D~ B^ —>• £)*+ Dj, B^ —>• £)*+ D*~ , or B^ — >• Df KJ, with clear signals, where I and vi 

stand for lepton and neutrino accordingly. The possibility of the production of the |(6s)[u]) or quarkonium 

via t-quark or t-quark decays is worth serious consideration especially when the LHC upgrades . 


C. Uncertainty analysis 

In the subsection, we discuss the uncertainties for the |(&s)[n.]) quarkonium production through top decays. For 
the present calculation, their main uncertainty sources include the nonperturbative bound state matrix elements, the 
CKM matrix elements, the renormalization scale (Xr, and the constituent quark masses rrib and mg. These parameters 
are the main uncertainty source for estimating the |(&s)[n]) quarkonium production. Here we only discuss the decay 
widths of the |(&s)[n]) quarkonium production through top quark decays under the constituent quark masses of the 
|( 6 s)[n]) quarkonium and the wave functions at the origin of the corresponding constituent quark masses in detail. In 
the following, we shall concentrate our attention on the uncertainties caused by iris and mb, whose center values are 
taken as mg = 0.50 ± 0.10 GeV and mb = 4.87 ± 0.20 GeV. And for clarity, when discussing the uncertainty caused 
by one parameter, the other parameters are fixed to be their center values. 

Typical uncertainties for mg, mb, and the values of the wave functions at the origin of the corresponding constituent 
quark masses are presented in Tables. 13 and ED In the Tables. El and ED it can be found that sizable uncertainties 
for varying mg and mb- The decay width will decrease with the increment mass of mg. But the decay width will 
increase with the increment mass of mb- And such the decay width will slow vary a tendency with a heavier quark 
mass. 

The uncertainties are drawn by shaded bands in the Figs. [S] and [T] In the Figs. [5] and [71 the up shaded band is the 
uncertainties of the |( 6 s)[l^S'i]) meson for the varying constituent quark masses, the middle shaded band is that of 
the |( 6 s)[l^S'o]) meson, the down shaded band is that of the |( 6 s)[lP]). In the above two shaded bands in the Figs. [5] 
and [71 where the solid-line in the center of the shaded bands is for mg = 0.50 GeV and mb = 4.87 GeV, the upper 
edge of the band is for mg = 0.40 GeV and mb = 4.67 GeV, which the correspondingly radial wave function at the 
origin is |7?|(bs)[is]) P = 0.4133 GeV^. And the lower edge of the bands is for mg = 0.60 GeV and mb = 5.07 GeV, 
and the correspondingly |i?|(bs)[is])P = 0.9404 GeV^. In the lowest shaded bands, the solid-line in the center of the 
shaded bands is for mg = 0.69 GeV and mb = 5.14 GeV of the |( 6 s)[lP]) constituent quark masses, the upper edge 
of the band is for mg = 0.59 GeV and mb = 4.94 GeV, and the correspondingly |i?|j^pj^(0)p = 1.999 x 10“^ GeV^, 
and the lower edge of the band is for mg = 0.79 GeV and mb = 5.34 GeV, and |i?|j^pj^(0)p = 4.472 x 10“^ GeV^, 
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TABLE V: Uncertainties for the decay width of the processes t —>■ |(fos)[n]) + W'^s under the B.T. potential (n/ = 3). 


S states 

ms (GeV) 
|P|( 6 .)[is]>|" iGeV^) 
r(t^ |(&s)[GSo])) (MeV) 
r(t^ |(6s)[l3S'i]))(MeV) 

0.40 

0.4162 

5.677 

11.41 

0.50 

0.6516 

4.528 

8.614 

0.6 

0.9316 

3.730 

6.463 


ms (GeV) 

0.59 

0.69 

0.79 


l^f(6.)[iP]>l" ( 10 -" GeU®) 

2.022 

3.028 

4.399 

P states 

r(t^|( 6 s)[liPi]))(KeV) 

433.6 

293.3 

214.7 


r(t^|( 6 s)[PPo]))(KeV) 

136.6 

119.7 

109.4 


r(t^|(6s)[l3Pi]))(KeV) 

203.8 

142.5 

107.6 


r(t^|( 6 s)[PP 2 ]))(KeV) 

363.1 

238.1 

169.0 


TABLE VI: Uncertainties for the decay width of the processes t —^ |(tis)[n]} + s under the B.T. potential (n/ =3). 


S states 

mb iGeV) 
|Pl( 6 .)[isi>l" (GeU^) 
Tit^\ibs)[l^So]}) (MeV) 
Tit^\ibs)[l^Si])) (MeV) 

4.666 

0.6458 

4.498 

8.468 

4.866 

0.6516 

4.528 

8.614 

5.066 

0.6569 

4.555 

8.751 


mb iGeV) 

4.94 

5.14 

5.34 


I^K 6 .)[iPi)l" ( 10 -" GeU^) 

2.985 

3.028 

3.070 

P states 

r(t^|(bs)[liPi]))(KeV) 

290.0 

293.3 

296.2 


r(t^|( 6 s)[PPo])) (KeV) 

125.6 

119.7 

114.2 


r(t^|( 6 s)[PPi])) (KeV) 

141.6 

142.5 

143.5 


r(t^|(te)[PP 2 ]}) (KeV) 

233.3 

238.1 

242.7 


where the contributions of |(6s)[lP]) form the color-singlet of |(6s)[l^Pi]) and |(6s)[l^Pj]) (with J = [0,1,2]) states 
have been summed up. 

Adding all the uncertainties caused by the constituent quark masses Wg = 0.50 ± 0.10 GeV and rub = 4.87 ± 0.20 
GeV in quadrature and the corresponding wave functions at the origin for the process t —>■ |(6s)[n]) -I- sW~^, we can 
obtain 


Tit ^ |( 6 s)[li,So]) + W+s) = 4528t“^®o KeV, 
Tit ^ |(6s)[1^5i]) + VU+s) = 8614^2799 KeV, 



FIG. 6 : (color online). Uncertainties of dV/ds\ and dV/ds 2 for t —^ |( 6 s)[n]) -|- sW^, where contributions from the color-singlet 
P-wave states have been summed up. The up shaded band, the middle shaded band, the down shaded band are for the 
uncertainties of |(fes)[l^S'i]), |(fes)[l^S'o]), and |( 6 s)[lP]), respectively. 
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FIG. 7: (color online). Uncertainties of dV/dcos 6\2 and dV/dcosOn for t —>■ |(&s)[n]) + sW^, where contributions from the 
color-singlet P-wctve states have been summed up. The up shaded band, the middle shaded band, the down shaded band are 
for the uncertainties of |(bs)[l®Si]), |(fos)[l^So]) and |(fos)[lP]}, respectively. 


|(6s)[lP])-bVU+s) = 793.6;?^^jKeV. (27) 

If the excited |(6s)[n]) quarkonium states decay to the ground spin-singlet S'-wave state |(6s)[l^S'o]) with 100% 
efficiency via electromagnetic or hadronic interactions, we can obtain the total decay width of the top quark decay 
channels under the B.T. potential. 

r{t^\{bs)[l^So]) + W+s) = 14.19t^;3® MeV. (28) 

IV. SUMMARY 

In the present paper, for studying the |(&s)[n.]) quarkonium production through t-quark or <-quark decays, we 
calculate the masses of the |(&s)[n]) quarkonium under the B.T. potential and the values of the Schrodinger radial 
wave function at the origin of the |(&s)[n]) quarkonium within the five potential models, and made a detailed study on 
the |(6s)[n]) quarkonium production via top quark and antitop quark semiexclusive decays channels, t —>■ |(6s)[n]) -I- 
W^s and i —|(&s)[?T']) + W~b, within the CSQCDFF. Results for the |(&s)[n]) quarkonium states, i.e., |(6s)[1^5o]), 
|(6s)[1^5'i]), |(6s)[l^Pi]), and |(6s)[l^Pj]) (with J = [0,1,2]) have been presented. And to provide the analytical 
expressions as simply as possible, we have adopted the “improved trace technology” to derive Lorentz-invariant 
expressions for top quark and antitop quark semiexclusive decay processes at the amplitude level. Such a calculation 
technology shall be very helpful for dealing with processes with massive spinors. 

Numerical results show that excited |(5s)[lP]) or |(6s)[lP]) states in addition to the ground IS'-wave states can 
also provide sizable contributions to the |(6s)[n]) or |(&s)[n]) quarkonium production through top quark decays, so 
one needs to take the excited states into consideration for a sound estimation. If all the excited states decay to the 
ground state |(6s)[1^5'o]), we can obtain the total decay width for the |(&s)) quarkonium production through top quark 
decays as shown by Eq. (PSI) . At the LHC, due to its high collision energy and high luminosity, sizable |(6s)[n]) or 
|(6s)[n]) quarkonium events can be produced in t-quark or t-quark decays; i.e., about 1.3 x 10® |(&s)) or |s[n])|(6s)) 
quarkonium events per year can be obtained. 

Acknowledgements: We are grateful to Professor Xing-Gang Wu for many enlightening discussions. This work 
was supported in part by the Natural Science Foundation of China under Grant No. 11347024, the Scientific and 
Technological Research Program of Chongqing Municipal Education Commission under Grant No. KJ1401313 and 
the Research Foundation of Chongqing University of Science and Technology under Grant No. CK2016Z03, the 
Natural Science Foundation Project of CQGSTG under Grant No. 2014jcyjA00030. 


V. APPENDIX 

A. The “improved trace technology” at the amplitude level 


We present its main idea of the “improved trace technology” at the amplitnde level. 
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First, we introduce a masses spinor with negative helicity M_(fco), which satisfies the following projection: 

U-{kQ)u-{kQ) = (29) 

where ko is an arbitrary light-like momentum, fcg = 0, ui- = (1 — 75 )/ 2 . Then, we construct the positive helicity 
states 


u+{ko) = ^iU-{ko). (30) 

where ki is an arbitrary space-like momentum, ki is defined as kf = —1 and ko ■ ki = 0 . u+{ko) is satisfied with the 
projection: 


u+{ko)u+{ko) = uj+^Q. 


(31) 


where a;+ = (1 -I- 75 )/ 2 . 

Using these two masses spinors, one can construct the massive spionrs for the fermion and antifermion as follows: 

Usiq) = (^ -I- m)u-{ko)/y/2ko ■ q, (32) 

u-s{q) = (^-l-m)u+(A:o)/\/2fc(I~g. (33) 


here the spin vector = ^ - -^kQ^. 

Using the above identities, one can write down the amplitude M±s±s' with four possible spin projections in the 
trace form for 7 -matrices: 

Mss> = 7Vrr[(|^-I-mt)a;_^Qcdot((^^-I-m^) • A)] 

M^ss' = NTr[{^ + mt ) uj +^ Q {{^^+ ms )- A)] 

M_ss' = NTr[{^ + ■ A)] 

Mss' = NTr[{$ + rnt)w+^^$Q{{(^^+'ms) ■ A)]. (34) 


where A = Yl, is the amplitudes of the above-mentioned decays channels. N = l/-\/4(fco • fc)(fco • qi) is the 

n—1 

normalization constant. 

It is easy check that M±s±8' is orthogonal for each other. Then, the squared amplitude M can be written as 


|M|2 = \Mss'? + \Ms-s'? + \Mss'? + \Mss'? 
= |Mi|2 + |M2|^ + |M3|2 + |M4|^ 


(35) 


where the four amplitude Mi with (i = 1, • • •, 4) are defined as 


Ml 

M 3 


Mss' + M-s-s' 

72 

Mss' - Mss' 

72 


,M 2 


,M4 


Mss' - Ms-s' 

72 

Ms-s' + M-ss' 

71 


(36) 


To simply the expression as much as possible, we take kg = qi — ak with a 
kuq 2 pqicr ( Nq ensures ki ■ ki = — 1 ), which leads to 


k-qi+y/ {k-qi)^-k^qf 

-P- 


and fci = 


h 


N 0 J 5 


k ■ 92^^ + <?1 • <? 2 ^ 


qi- kqfe- 


Then, the resultant Mi are. 

Ml = LixTr[{ll: + mt){^^+ms)A] 

M2 = -L2'KTr[{$ + mtYs{i^-^ms)A\ 

M 3 = M 3 / -I- iVo[(/c • 93)7715 - (qi • (?3)mt]M2 
M 4 = M 4 < - iVo[(/c • 93 ) 7775 -I -(91 • 93 )? 77 t]Mi 


(37) 

(38) 

(39) 

(40) 
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where L 1.2 = l/(2i/fc • qi ± mtiris) and 



—Tr + mt)(^{^^+ms)A 


Tr + mt)"f 5 (^i^^+ms)A 


(41) 


(42) 
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